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Abstract 



We report the results of a search for point-like deviations from isotropy in the 
arrival directions of ultra-high energy cosmic rays in the northern hemisphere. In 
the monocular data set collected by the High-Resolution Fly's Eye, consisting of 
1,525 events with energy exceeding 10^^'^ eV, we find no evidence for point-like 
excesses. We place 90% c.l. upper limits less than or equal to 0.8 cosmic rays/km^yr 
on the flux from such sources as a function of position in the sky. 
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1 Introduction 

In the search for the origins of ultra-high energy cosmic rays (UHECR), a 
variety of strategies have been employed to detect deviations from isotropy in 
event arrival directions. The results of these searches have been ambiguous. 
The most recent results [1] have supported the null hypothesis for large-scale 
dipolc behavior in arrival directions for particles above 10^^'^ eV. However, 
previous experiments [2] have found evidence for a dipole moment in Right 
Ascension (RA) at energies above 10^^ eV. On smaller angular scales, ex- 
cesses have been alternately claimed and refuted in the vicinity of Cygnus X- 
3 [3,4,5,6], an x-ray binary within our galaxy, including the report of a possible 
excess in a point-source search [7] . 

Point-like excesses at these energies can arise from only a limited number of 
source scenarios. Galactic and extragalactic magnetic fields are expected to 
produce large perturbations in the arrival directions of charged particles: A 
proton of energy 10^^'^ eV may be deflected by several tens of degrees as it 
traverses the disk of the Milky Way galaxy, with a typical magnetic field of 
order 1 microgauss [8]. A compact arrival direction excess at these energies 
would therefore suggest neutral primaries. Neutrons however possess a fifetime 
of 3 X 10^^ seconds at 10^^'^ eV, and therefore cannot have originated more 
than 30 kpc from Earth. Thus any viable source of standard model neutral 
hadronic matter would have to be located within the Milky Way Galaxy. 

In this paper, we conduct a search for point-like behavior in arrival direction of 
cosmic ray events above 10^^'^ eV in the northern hemisphere, using a skymap 
technique in which we evaluate our sensitivity using Monte Carlo simulated 
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sources. In addition, we consider the historically significant source candidate 
Cygnus X-3 as the focus of an a priori search. 



2 The HiRes Monocular Data Set 

The High-Resolution Fly's Eye (HiRcs) consists of two nitrogen fluorescence 
observatories — HiRes-1 and HiRes-11 — separated by 12.6 km and located 
at Dugway, Utah. HiRes was conceived as a stereo detector, however due to 
the larger available statistics it is desirable to reconstruct extensive airshow- 
ers in monocular mode as well. This HiRes-1 monocular data set consists of 
2,820 good-weather detector hours of data, collected between May 1997 and 
February 2003. A total of 1,525 events with energies exceeding 10^^'^ eV were 
collected during this time and are included in the present analysis. 

The HiRes-1 monocular data set and airshower reconstruction by profile con- 
strained fitting has recently been described in the literature [9]. A residual 
effect of the profile-constrainted fit technique is orientation-dependent (ellip- 
tical) uncertainties in the airshower arrival directions. In Figure 1, the air- 
shower reconstruction geometry is illustrated for a monocular air fluorescence 
detector. In this view, the shower-detector plane (SDP) for HiRes-I events is 
well-reconstructed, with uncertainty parameterized as 

(TSDP = 88.151°e-^^^^ + 0.374° (1) 

where Ax is the angular tracklength of the shower in degrees. Typical values 
of (7 SDP for this analysis range from 0.4° — > 1.7°. The angle of the track within 
the SDP, is less well reconstructed and is parameterized by 

(7^ = 18.381°e- i°gio(i-45E) ^ 4 073° (2) 

where the energy E is expressed in EeV (lO-*^^ eV). Typical values of cr* in this 
analysis range from 5.4° 15°. 

In Figure 2, we plot the skymap formed from the arrival directions of events in 
the present data set. Each event's "error ellipse" is represented by generating 
1,000 points per event, distributed according to the Gaussian error model of 
Equations 1 and 2. Figure 2 is plotted in equatorial coordinates as a polar 
plot. Note that bins are assigned using a cartesian projection of the polar 
plot shown in Figure 2 and all similar figures. As such, angular bin size varies 
across the map, but averages approximately 1° x 1°. 

We next discuss the Monte Carlo technique by which we evaluate the signif- 
icance of fiuctuations in the skymap as well as our sensitivity to point-like 
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Fig. 1. The geometry of reconstruction for a monocular air fluorescence detector, 
behavior in arrival direction. 
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Fig. 2. Skymap of arrival directions of events in the the HiRes-I monocular data set, 
plotted in polar projection, equatorial coordinates. Each HiRes event is represented 
by 1,000 points randomly thrown according to the elliptical Gaussian error model of 
Equations 1 and 2. The bin size in this plot (and all similar plots) is approximately 
1° X 1°. 
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3 The Monte Carlo; Compcirison of Data to Expectation from an 
Isotropic Background 



We use a library of simulated events, generated by the Monte Carlo tech- 
nique and reconstructed using the profile-constrained reconstruction program 
to determine the background expectation for isotropically distributed sources 
as well as to evaluate our sensitivity to point-like behavior in arrival direction. 
For this library, an isotropic distribution is assumed for events possessing the 
spectrum and composition reported by the stereo Fly's Eye experiment [10,11]. 

A detector runtime database is used to randomly assign a time from detector 
"on" periods to each event in the isotropic background data set. A total of 
1,000 isotropic data sets with the same sky exposure as the HiRes-I monocular 
data set were generated for comparison studies. Further discussion of this 
Monte Carlo can be found in Reference [1]. In Figures 3 and 4 we compare the 
data and Monte Carlo distributions of events in the variables RA and DEC, 
respectively. 
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Fig. 3. Comparison of HiRes-I data (points) and Monte Carlo (solid histogram) 
distributions in right ascension (RA). 

In order to understand the significance of the fiuctuations in Figure 2, we 
compare the data on a bin-by-bin basis to the 1,000 simulated data sets. 
Defining No at a as the bin density of the data, Nmc as the bin density of the 
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Fig. 4. Comparison of HiRes-I data (points) and Monte Carlo (solid histogram) 
distributions in declination (DEC). 

simulated isotropic data sets, and auc ^ the standard deviation of the Monte 
Carlo bin density, the variable 



Xi 



{Ndata - (Nmc)) 



(3) 



provides a measure of the fluctuation per bin. Figure 5 shows the distribution 
of Xi as a function of position in the sky for the HiRes-I monocular data set 
as extracted from this technique. 

The bin-by-bin distributions of Xi ^re non-Gaussian (Figure 6) and vary as 
a function of position in the sky. Thus it is necessary to develop a technique 
to evaluate the significance of possible sources. Our technique uses the Xi in- 
formation in neighboring bins to pick out significant fiuctuations above back- 
ground from the skymap. The parameters in the technique are tuned on sim- 
ulated point-like sources. 
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Fig. 6. Example distribution of xi (Equation 3) values for 1,000 MC data sets in 
the bin located at 5 hours RA, 40° DEC. 
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4 The Monte Carlo; Simulation of Point-Like Sources 

We have two objectives in simulating point-like sources: The first is using these 
simulated sources to tunc point source selection criteria. Secondly, simulated 
sources provide a straighforward method by which to quantify our sensitivity 
to point-sources and derive flux upper limits. 

Simulated source skymaps are created by randomly replacing events in a 
simulated isotropic data set with A^^ events at the chosen position for the 
source. The central-value coordinates of the simulated source event are ran- 
domly shifted according to the error ellipse, which is taken from the replaced 
event. The shift simulates the effect of detector resolution only. Finally, the 
orientations of error ellipses are randomized. 

An example of a simulated source is shown in Figure 7. This source is superim- 
posed on a Monte Carlo data set in Figure 8, and the quantity xi (Equation 3) 
is evaluated for each bin in Figure 9. We note that source events overlap in 
a fairly small angular region, as seen in Figure 7. Thus we have sensitivity to 
fairly compact deviations from isotropy, in spite of our elongated error ellipses. 
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Fig. 7. An Ns = 25 event source shown without the isotropic background. The 
source has been inserted at 5 hours RA, 40° DEC. Each source event is represented 
by 1,000 points randomly thrown according to the elhptical Gaussian error model 
of Equations 1 and 2, where the error ellipse is taken from the replaced isotropic 
event. 
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Fig. 8. Skymap of arrival directions of events for a simulated data set, having the 
same overall exposure as the HiRes-I monocular data set, with a 25 event source 
superimposed at 5 hours RA, 40° DEC (compare to Figure 2). 
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5 Calculation of Significances 



We now describe a procedure by which we can identify point-hke behavior in 
arrival direction (for example, the simulated source of Figures 7, 8, and 9) 
while simultaneously rejecting false positives arising from fluctuations of the 
background. 

Due to detector resolution, it is desirable that we search for sources by consid- 
ering points over an extended angular region. We consider a "search circle" of 
radius R, where R is expressed as an angle in degrees. Within the search circle, 
we count the fraction of bins F having a Xi value greater than some threshold 
Xthr- The parameters R and Xthr are chosen to optimize the signal size, and 
a cut is chosen on the fraction F which reduces the false positive probability 
to an acceptable level. 

Our maximum sensitivity to point-like behavior in arrival direction, given the 
HiRes-I pointing uncertainty, was determined to require a search circle of 
R — 2.5°, and a value Xthr — 4. (In the case in which the bin densities 
are normally distributed, this corresponds to 4a.) The optimum values for 
these parameters were determined by simulating sources at various locations 
in the sky and maximizing our sensitivity to these sources. The values for these 
parameters are found to be largely insensitive to the position in the sky and 
the number of events in the source. Additionally, small variations in either of 
these parameters do not have a significant impact on our results. 

Due to low statistics at the edge of HiRes' acceptance, we consider only search 
circles with centers whose declinations are greater than 0°. That is, we only 
search for sources north of the celestial equator. Approximately 10% of HiRes 
events have central-value coordinates south of the equator. These events can 
contribute to the search if their error ellipses extend north of DEC = —2.5°. 

In Figure 10, we have plotted for each bin the fraction F, for R — 2.5° and 
Xthr = 4, of the simulated point source of Figure 8. The simulated source 
stands out clearly in this figure. 

The final parameter in this search algorithm is the cut placed on the quan- 
tity F. We evaluate this cut by requiring that the probability of a simulated 
isotropic data set - without a superimposed simulated source — exceeding the 
cut be no more than 10% over the entire sky (Figure 11). We choose a cut 
value oi F — 0.33, corresponding to a false-positive probability of 10%. 

Figure 12 shows the F distribution for the monocular data set. The "hottest" 
spot on this graph, near DEC = 20° and RA = 20 hours, has a value F = 
0.15. 87% of simulated isotropic data sets have F > 0.15 (see Figure 11). We 
conclude that our observation is consistent with a fluctuation from an isotropic 
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Fig. 10. F distribution derived from the xi ™ap of Figure 9. F is the fraction of 
bins within radius R = 2.5° having a xi value of 4 or greater. The simulated source 
at 5 hours RA, 40° DEC clearly stands out as having an exceptionally large excess 
fraction F. 

background. 



Next, we evaluate the corresponding sensitivity and flux upper limits as a 
function of position in the sky. 
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Fig. 11. Occurrence rate of false positives versus F, for a 2.5° search circle and xi 
threshold of 4. A cut at F = 0.33 corresponds to a false-positive probability of 10%. 
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Fig. 12. F distribution derived from the HiRes-I monocular data (Figures 2 and 
5). i*" is the fraction of bins within radius R = 2.5° having a xi value of 4 or greater. 
No points satisfy our search criteria. The largest F occurs at DEC = 20° and RA 
= 20 hours and has a value F = 0.15, corresponding to a false-positive probability 
of 87%. 
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6 Sensitivity and Upper Limits 



In the preceding section, we found no evidence for the presence of point-hke 
excesses in the HiRes-I monocular datasct above 10^^'^ eV. Further, the signifi- 
cance of the largest point-like fluctuation in the data is well below the threshold 
established to minimize the likelihood of false positives in an isotropic distri- 
bution. To quantify our null result, we follow the suggestion of Feldman and 
Cousins [12] and calculate both a set of flux upper limits and the "sensitivity" 
of the experiment to such point-like excesses. The results of these calculations 
are reported below. 

6. 1 Sensitivity 

The sensitivity of the experiment is defined as the average 90% confidence level 
flux upper limit that would be reported by an ensemble of like experiments 
with no true signal. Since this average upper limit will vary as a function of 
position on the skymap due to different background expectations, we calculate 
our sensitivity at set of gridpoints (Table 1) distributed evenly across the 
Northern Hemisphere. We choose the right ascension values of our gridpoints 
to correspond approximately to the HiRes "solstices" and "equinoxes" , i. e. to 
the RA lines of high/low and midrange event statistics. 

To determine our sensitivity to a number of "source" events (A''^) at a given 
gridpoint, we generate 400 simulated isotropic datasets with point-like sources 
superimposed. The number of events in each source is Poisson distributed with 
mean value {Ns). These datasets contain the same total number of events as 
the HiRes-I data. 

We then determine the percentage of trials at each location for which our 
reconstruction algorithm "finds" a source of size {Ng). In the case of the 
sensitivity calculation, wc say the algorithm "finds" a source if at one point 
on the skymap F fluctuates above our preselected threshold value of F = 0.33. 
The value of {Ns) for which signal was declared for 90% or better of the trials 
was termed 

The distribution of A^33 at our grid points is illustrated in Figure 13. Wc 
estimate that the systematic uncertainty in the calculation of A 33, due to 
uncertainties in the size of the error ellipses, is < 1 event. 

The HiRes-I detector flux sensitivity at each grid point is A^.33 at that point di- 
vided by the local exposure. We calculate the detector exposure [13] for point 
sources at the grid points by the following procedure: Monte Carlo events are 
generated at the grid points, assigned a time from the distribution of HiRes 
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Fig. 13. Numerical values of — the mean number of source events for which 
signal was declared in 90% or better of 400 trials with a cut on F of .33 — at 21 
grid points in the Northern Hemisphere. Exact numbers and locations of gridpoints 
are given in Table 1. The systematic uncertainty in the calculation of A^.33, due to 
uncertainties in the size of the error ellipses, is < 1 event. 



detector ontimes, and projected towards the detector aperture. Local coordi- 
nates and times are determined, then the event is paired with a shower from 
the Monte Carlo event library having similar local coordinates. An attempt is 
then made to reconstruct the Monte Carlo event with the profile-constrained 
fitting technique. The exposure, defined as the fraction of events reconstructed 
multiplied by the detector aperture (area) and time, can then be used to deter- 
mine flux sensitivity (as well as flux upper limits) for each of the grid locations. 
These exposures are listed in Table 1. The final fiux sensitivities are shown in 
Table 1 and Figure 14. 
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Fig. 14. Flux sensitivity (cosmic rays/km^yr) at 21 grid points in the Northern 
Hemisphere. Sensitivity is calculated by dividing (Figure 13) by the exposure 
at each grid point. The numbers along with the exact locations of grid points are 
given in Table 1. 

6.2 Upper Limits 



We place our 90% confidence level flux upper limits by making use of the 
fact that F never fluctuated above 0.15 in the HiRes-I data. (Figure 10. We 
determine the value of — the mean number of source events for which 
signal was declared in 90% or better of trials with a threshold of F = .15 — 
in the same manner in which we determined A^33 in the previous section, he 
results are summarized in Figure 15. Our 90% c.l. flux upper limit at each 
grid point (Figure 16) is at that point divided by the local exposure. The 
largest flux upper limit across the entire sky is 0.8 cosmic rays above 10^^-^eV 
per km^yr. 

Finally, we note that the a priori source candidate in the direction of Cygnus 
X-3 (RA 20.5 hours, DEC 40.7°) is very near the grid point located at RA 
20.5 hours, DEC 45°. We place a 90% c.l. flux upper limit from Cygnus X-3 at 
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0.5 cosmic rays above lO-'^^'^eV per km^yr. Previous Cygnus X-3 flux results 
were drawn from events samples with energies above 5 x 10^^ cV, so a direct 
comparison is impossible. We can infer — assuming that any cosmic rays from 
Cygnus X-3 have an energy spectrum similar to that of the full sky — that 
an extrapolation of our result is not competitive with prior upper limits for 
events above 5 x 10^^ eV. However, this is the first reported measurement for 
a high-statistics sample above 10^^'^eV. 




Fig. 15. Numerical values of iV.15 — the mean number of source events for which 
signal was declared in 90% or better of 400 trials with a cut on F = .15 — at 21 
grid points in the Northern Hemisphere. Exact numbers and locations of gridpoints 
are given in Table l.The systematic uncertainty in the calculation of iV.15, due to 
uncertainties in the size of the error ellipses, is < 1 event. 



21 




Fig. 16. 90% c.l. Flux Upper Limit (cosmic rays/km^yr) at 21 grid points in the 
Northern Hemisphere. Flux upper limit is calculated by dividing by the ex- 
posure at each grid point. Exact numbers and locations of gridpoints are given in 
Table 1. 
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Table 1 

Locations of gridpoints, threshold signal strengths and A^.15, exposures (with 
uncertainty 5%, primarily from Monte Carlo statistics), detector flux sensitivity, 
and 90% confidence level flux upper limits, for cosmic rays with energy exceeding 
10^^-5 eV. 
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7 Conclusions 



We have conducted a search for point-hke excesses in the arrival direction of 
ultra-high energy cosmic rays with energy exceeding 10^^'^ cV in the northern 
hemisphere. We place an upper limit of 0.8 cosmic rays/(km^ yr) (90% c.l.) 
on the flux from such sources across the entire sky and place more stringent 
limits as a function of position. We also determine sensitivity as a function of 
position in the sky. The HiRes-I monocular data is thus consistent with the 
null hypothesis for point-like excesses of neutral primary cosmic rays in this 
energy range. 
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